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Abstract 

The reductive perturbation method has been employed to derive 
the Korteweg-de Vries (KdV) equation for small but finite amplitude 
electrostatic ion-acoustic waves in unmagnitized collisionless weakly 
relativistic warm plasma. The Lagrangian of the time fractional KdV 
equation is used in similar form to the Lagrangian of the regular KdV 
equation. The variation of the functional of this Lagrangian leads to 
the Euler-Lagrange equation that leads to the time fractional KdV 
equation. The Riemann-Liouvulle definition of the fractional deriva- 
tive is used to describe the time fractional operator in the fractional 
KdV equation. The variational-iteration method given by He is used 
to solve the derived time fractional KdV equation. The calculations 
of the solution with initial condition Aq sech(cx)^ are carried out. The 
result of the present investigation may be applicable to some plasma 
environments, such as ionosphere. 

Keywords : Ion-acoustic waves; Euler-Lagrange equation, Riemann- 
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1 Introduction 

Nonlinear evolution equations are widely used as models to describe com- 
plex physical phenomena. Because most classical processes observed in the 
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physical world are nonconservative, it is important to be able to apply the 
power of variational methods to such cases. A method used a Lagrangian 
that leads to an Euler-Lagrange equation that is, in some sense, equivalent 
to the desired equation of motion. Hamilton's equations arc derived from 
the Lagrangian and are equivalent to the Euler-Lagrange equation. If a La- 
grangian is constructed using noninteger-order derivatives, then the resulting 
equation of motion can be nonconservative. It was shown that such fractional 
derivatives in the Lagrangian describe nonconservative forces [1, 2]. Further 
study of the fractional Euler-Lagrange can be found in the work of Agrawal 
[3, 4], Baleanu and coworkers [5, 6] and Tarasov and Zaslavsky [7, 8]. During 
the last decades, Fractional Calculus has been apphed to almost every field 
of science, engineering and mathematics. Some of the areas where Fractional 
Calculus has been applied include viscoclasticity and rhcology, electrical engi- 
neering, electrochemistry, biology, biophysics and bioengineering, signal and 
image processing, mechanics, mechatronics, physics, and control theory [9]. 

The propagation of solitary waves is important as it describes charac- 
teristic nature of the interaction of the waves and the plasmas. In the case 
where the velocity of particles is much smaller than that of light, ion-acoustic 
waves present the non-relativistic behaviors, but in the case where the ve- 
locity of particles approaches that of light, the relativistic effect becomes 
dominant [10]. Actually high speed and energetic streaming ions with the 
energy from 0.1 to 100 MeV are frequently observed in solar atmosphere 
and interplanetary space. Nevertheless, relativistic ion-acoustic waves have 
not been well investigated. When we assume that the ion energy depends 
only on the kinetic energy, such plasma ions have to attain very high veloc- 
ity of relativistic order. Thus, by considering the weakly relativistic effect 
where the ion velocity is about ^ of the velocity of hght, we can describe the 
relativistic motion of such ions in the study of nonlinear interaction of the 
waves and the plasmas [11]. It appears that the weakly relativistic and ion 
temperature effects play an important role to energetic ion-acoustic waves 
propagating in interplanetary space [12]. Washimi and Taniti [13] were the 
first to use reductive perturbation method to study the propagation of a slow 
modulation of a quasimonochromatic waves through plasma. And then the 
attention has been focused by many authors [14-15]. 

To the author's knowledge, the problem of time fractional KdV equation 
in weakly relativistic plasma has not been addressed in the literature before. 
So, our motive here is to study the effects of time fractional parameter on 
the electrostatic structures for a system of collisionless plasma consisting of a 
mixture of warm ion-fiuid and isothermal electrons with ion flow velocity has 
a weakly relativistic effect. We expect that the inclusion of time fractional 
parameter and a weakly relativistic effect will change the properties as well 



2 



as the regime of existence of solitons. 

Several methods have been used to solve fractional differential equations 
such as: the Laplace transform method, the Fourier transform method, the 
iteration method and the operational method [16]. Recently, there are some 
papers deal with the existence and multiplicity of solution of nonlinear frac- 
tional differential equation by the use of techniques of nonlinear analysis 
[17-18]. In this paper, the resultant fractional KdV equation will be solved 
using a variational-iteration method (VIM) firstly used by He [19]. 

This paper is organized as follows: Section 2 is devoted to describe the for- 
mulation of the time-fractional KdV (FKdV) equation using the variational 
Euler-Lagrange method. In section 3, Variational-iteration Method (VIM) is 
discussed. The resultant time-FKdV equation is solved approximately using 
VIM. Section 5 contains the results of calculations and discussion of these 
results. 



2 Basic equations and KdV equation 

Consider coUisionless ionization-free unmagnetized plasma consisting of a 
mixture of warm ions-fluid and isothermal electrons. Assume that the ions 
flow velocity has a weak relativistic effect, and therefore there exist stream- 
ing ions in an equilibrium state when sufficiently small- but finite-amplitude 
waves propagate one-dimensionally. Such a system is governed by the fol- 
lowing normalized equations [12]: 

d d 

— n(x, t) + t)u{x, t)] = 0,(la) 

[| + .(x,^)^][7(x,^)K(x,^)] + ^|p(x,^) + ^0(x,^) = 0,(lb) 

[^ + u{x,t)^]p{x,t)+3p{x,t)^[^u{x,t)] = 0,(lc) 

—4>{x,t)+n{x,t)-ne{x,t) = 0.(ld) 

The electrons temperature Tg is much larger than the ions temperature 
Tj and in this case, for simplicity, one can neglect the inertia of the electrons 
relative to that of the ions, i.e. the high-frequency plasma oscillations are 
neglected. Since we are interested in the regime of density and velocity 
fluctuations near the ion plasma frequency, so the isothermal electrons given 
by: 

ne{x,t) = cxp[(l){x,t)], (le) 
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for weakly relativistic effects, the relativistic factor 7(2;, t) — [l—u{x, tY / Cq] 2 
is approximated by 

7(x,t)^l + tz(x,t)V(2c2), (If) 

where c is the velocity of light. In equations (1) n{x,t) and ne{x,t) are 
the densities of ions and electrons, respectively, u{x,t) is the ions flow ve- 
locity, p{x, t) is the ions pressure, (f){x, t) is the electric potential, x is the 
space co-ordinate and t is the time variable, a — Ti/T^. « I is the ratio 
of the ions temperature to the electrons temperature. All these quantities 
are dimensionless, being normalized in terms of the following characteristic 
quantities: n{x,t) and ne{x,t) are normalized by the unperturbed electrons 
density no, u{x,t) and c are normalized by the sound velocity {kBT^./nii)^^'^ , 
p{x,t) and (j){x,t) are normalized by n^kBTi and kBT^/e, respectively and 
t and X are normalized by the inverse of the plasma frequency uj~l — 
(47re^no/mi)~^/^ and the electron Debye length Ad = (/csTe/47re^no)^^^, re- 
spectively, ks is the Boltzmann's constant, e is the elecron charge and is 
the mass of plasma ion. 

According to the general method of reductive perturbation theory, we 
introduce the stretched variables 

T^eh i^e^{x-Xt), (2) 

where A is the the phase velocity. All the physical quantities appeared in (1) 
are expanded as power series in terms of the amplitude of the perturbation 
e about the equilibrium values as 

n(e,r) = l + eni(e,r) + eV(e,T) + ..., (3a) 

u{i,T) = iio + eMi(C,r)-FeV(C,7') + --, (3b) 

p(e,r) = l + epi(e,r) + e2p2(e,T) + ..., (3c) 

0(e,T) = e(/)i(e,T) + e^02(e,T) + ..., (3d) 

with the boundary conditions that as|^| ^oo,n = ne=p=l,M = Mo, 

Substituting (2) and (3) into the system of equations (1) and equating 
coefficients of like powers of e. Then, from the lowest, second-order equa- 
tions in e and by eliminating the second order perturbed quantities n2(^, r), 
W2('C, t), P2(C) ''') and 02(C) ''')) we obtain the following KdV equation for the 
first-order perturbed potential: 

-0i(e, t) + A 0,(e, r)^0i(e, t) + B ^0,(e, r) = 0, (4a) 
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where 



^ = A-(A=-3.)2i + |(3T,-l), B = ^^^, (4b) 
with the transendental equation for A as 



l_(A2-3a)7i = =^ X^±J^ni±l^ (4c) 



and 



3 „o 3 R „ tin , 

71 = 1 + ;^^?', 12 = 7;- R=-- (4d) 

2 2 Co Co 

In equation (4a), (f>i{C, t) is a field variable, a space coordinate in 

the propagation direction of the field and r G T (= [0,To]) is the time.f 

The resultant KdV equation (4a) can be converted into time-fractional 
KdV equation as follows: Using a potential function V"(^, r) where 0i^(^, r) = 
^(C) ''") = ^(C) 7") gives the potential equation of the regular KdV equation 
(1) in the form 

V^ri^, t) + A V^{^, t)v^^{^, t)+B %^^(e, r) = 0, (5) 

where the subscripts denote the partial differentiation of the function with 
respect to the parameter. The Lagrangian of this regular KdV equation (4a) 
can be defined using the semi-inverse method [20, 21] as follows. 

The functional of the potential equation (5) can be represented by 

J{v)^ Jd^J dT{y(e,T)[ciy^,(e,r)+c2^y^(e,T)t;^e(e,T)+c3We^^^(e,T)]}, 

R T 

(6) 

where Ci, C2 and C3 are constants to be determined. Integrating by parts and 
taking = V^\r = V^\t = lead to 

J{v)^ Jd^J dT{y(e,T)[-ciye(e,T)K(e,T)-^c2Ay/(e,T)+c3W|(e,T)]}. 

R T 

(7) 

The unknown constants Cj (?' = 1, 2, 3) can be determined by taking 
the variation of the functional (7) to make it optimal. Taking the variation 
of this functional, integrating each term by parts and making the variation 
optimum give the following relation 

2ciV^ri^, r) + 3c2AV^{^, r)%(e, r) + 2csBV^^^^iC, r) = 0. (8) 
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As this equation must be equal to equation (5), the unknown constants 
are given as 

ci = 1/2, C2 = 1/3 and C3 = 1/2. (9) 

Therefore, the functional given by (7) gives the Lagrangian of the regular 
KdV equation as 

L(K, = -^^e(e,r)K(e,r) - + ^W|(e,r). (10) 

Similar to this form, the Lagrangian of the time-fractional version of the 
KdV equation can be written in the form 

< q;<1, (11) 

where the fractional derivative is represented, using the left Riemann-Liouville 
fractional derivative definition as [16] 

k-l < a<k,te[a,h]. (12) 
The functional of the time-FKdV equation can be represented in the form 

J{V) = j dij drFioD^V, y^,%), (13) 

R T 

where the time-fractional Lagrangian F{qD^V, V^, V^^) is defined by (11). 

Following Agrawal's method [3, 4], the variation of functional (13) with 
respect to V{^,t) leads to 

sj(v) = j,(j M^WV + l^^n + §-w,,Y (14) 

R T 

The formula for fractional integration by parts reads [3, 16] 
fdt fit) ,Dtg{t) = fdt g{t) tD^fit), fit), g{t) e [a, b], (15) 

J a J a 

where tD^, the right Riemann-Liouville fractional derivative, is defined by 
[16] 

k-l < a<l,te[a,b]. (16) 
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Integrating the right-hand side of (14) by parts using formula (15) leads 

to 

R T 

where it is assumed that 6V\t = SV\ii = SV^\r = 0. 

Optimizing this variation of the functional J{V), i. e;5J(V") = 0, gives 
the Euler- Lagrange equation for the time-FKdV equation in the form 

'^?»(^>-5fW^^W^"- '''' 

Substituting the Lagrangian of the time-FKdV equation (11) into this 
Euler-Lagrange formula (18) gives 

-\ r^?on(e,T) + ^ oD-,V^{i,T)+AV^{i,T)V^^{i,T)+BV^^^^{i,T) = 0. (19) 

Substituting for the potential function, V^{^,t) = 0i(C,t) ~ $(^,r), 
gives the time-FKdV equation for the state function $(^, r) in the form 

\[^D'^Mi, r) -r D^^, r)] + A $(e, r) ^^{^ r) + B %^(e, r) = 0, (20) 

where the fractional derivatives qD^ and tD^^ are, respectively the left and 
right Riemann-Liouville fractional derivatives and are defined by (12) and 
(16). 

The time-FKdV equation represented in (20) can be rewritten by the 
formula 

^ ^D^m t)+A $(e, r) $^(^, t) + B $^^^(e, r) = 0, (21) 

where the fractional operator qD'^ is called Riesz fractional derivative and 
can be represented by [4, 16] 

k-1 < a<l,t e[a,b]. (22) 

The nonlinear fractional differential equations have been solved using dif- 
ferent techniques [16-20]. In this paper, a variational-iteration method (VIM) 
[21] has been used to solve the time-FKdV equation that formulated using 
Euler-Lagrange variational technique. 
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3 Variational- iteration method 



Variational-iteration method (VIM) [21] has been used successfully to solve 
different types of integer nonlinear differential equations [22, 23]. Also, VIM 
is used to solve linear and nonlinear fractional differential equations [24, 25]. 
This VIM has been used in this paper to solve the formulated time-FKdV 
equation. 

A general Lagrange multiplier method is constructed to solve non-linear 
problems, which was first proposed to solve problems in quantum mechanics 
[21]. The VIM is a modification of this Lagrange multiplier method [22]. The 
basic features of the VIM are as follows. The solution of a hnear mathemat- 
ical problem or the initial (boundary) condition of the nonlinear problem is 
used as initial approximation or trail function. A more highly precise approx- 
imation can be obtained using iteration correction functional. Considering 

a nonlinear partial differential equation consists of a linear part LU{x,t), 

nonlinear part NU{x,t) and a free term f{x,t) represented as 

LU{x,t) + NU{x,t) = f{x,t), (23) 

where L is the linear operator and is the nonlinear operator. According 
to the VIM, the {n+ l)th approximation solution of (23) can be given by the 
iteration correction functional as [24, 25] 

Un+i{x,t)^Un{x,t)+ [ dTX{r)[LUn{x,r) + NUn{x,r)-f{x,r)], n>0, 

Jo 

(24) 

where A(r) is a Lagrangian multiplier and f/„(a;,r) is considered as a re- 
stricted variation function, i. e; 6Un{x,T) = 0. Extreme the variation of 
the correction functional (24) leads to the Lagrangian multiplier A(t). The 
initial iteration can be used as the solution of the hnear part of (23) or the 
initial value U{x,0). As n tends to infinity, the iteration leads to the exact 
solution of (23), i. e; 

U{x,t) = limUn{x,t). (25) 

n—^oo 

For linear problems, the exact solution can be given using this method 
in only one step where its Lagrangian multiplier can be exactly identified. 
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4 Time-fractional KdV equation solution 

The time-FKdV equation represented by (21) can be solved using the VIM 
by the iteration correction functional (24) as follows: 

Affecting from left by the fractional operator on (21) leads to 



dr ' u-r -v^''^i--r(a-l) 

8 8^ 
- ^D^'-^iA m r)g-^m r) + B r)], 

< «<1, re [0,To], (26) 
where the following fractional derivative property is used [16] 

{t - a) -"-J' 



k-l < P <k. (27) 



As a < 1, the Riesz fractional derivative qD^'^ ^ is considered as Riesz 
fractional integral that is defined by [16] 



1 1 



2Via) 



cir|t-r|"-7(r), a > 0. (28) 



where Ql^°'f{t) and rh"'fi^) right Riemann-LiouvuUe frac- 

tional integrals, respectively [16]. 

The iterative correction functional of equation (26) is given as 



$,+i(e,r) = $„(^,r) + 2'dT'A(r'){^$„(e,r') 



-^/^r"$„(e,r')|.'=o 



r 



rfa- 1) 



8 ~ 8^ ~ 

where n > and the function $„(^,r) is considered as a restricted variation 
function, i. e; 5$„(^, r) = 0. The extreme of the variation of (29) using the 
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restricted variation function leads to 



= r) + A(r) r) - £ c?r'^A(r') r') = 0. 

This relation leads to the stationary conditions 1+A(r) = and ^A(t') = 
0, which lead to the Lagrangian multiplier as A(r') = —1. Therefore, the 
correction functional (29) is given by the form 



'0 



T{a-1) 

+ ^DJ--[A $n(e,r')^3>n(e,r') + 5 ^*n(e,r')]}(30) 

where n > 0. 

In Physics, if r denotes the time-variable, the right Riemann-Liouville 
fractional derivative is interpreted as a future state of the process. For this 
reason, the right-derivative is usually neglected in applications, when the 
present state of the process does not depend on the results of the future 
development [3]. Therefore, the right-derivative is used equal to zero in the 
following calculations. 

The zero order correction of the solution can be taken as the initial value 
of the state variable, which is taken in this case as 

MC, r) = m 0) = ^0 sec h(cO', (31) 

where Aq = ^ and c = are constants and v is the travelling wave 

velocity. 

Substituting this zero order approximation into (30) and using the defi- 
nition of the fractional derivative (22) lead to the first order approximation 
as 

$i(e,r) = Aosech(ce)' 

+2Aoc sinh(cO sech(cO^ 

*[4c^B + {AoA - 12c'E)sech(cO'] ^, (32) 

L[a + 1) 



10 



Substituting this equation into (30), using the definition (22) and the 
Maple package lead to the second order approximation in the form 

MC,r) = Aosech{cCf 

+2Aoc sinh(cO sech(cO^ 

*[4c^B + {AqA - l2c^B) sec \i{c^f' ^ 



+2Aoc^ sech(cO^ 

*[32c^B^ + lQc^B{bAQA - 63c^B) sech(cO^ 
+2{3AlA^ - 176Aoc^AB + IGSOc^S^) sech(cO^ 

-7{AlA^ - A2Ao^AB + SeOc^S^) sech(ce)^]— -— - 

i (za + Ij 

+4AoC^ sinh(cO sec h(cO^ 
*[32c^52 + 2Ac^B{AqA - lAc^B) sech(cO^ 
+A{AlA^ - 32Aoc^AB + 2A0c^B^) scch(cO^ 
-5{AlA^ - 2AA^^c^AB + lAAc^B"^) sech(cOl 

r(2a + 1) r'^" 
*[r(a + l)]2r(3a + l)' 



(33) 



The higher order approximations can be calculated using the Maple or 
the Mathematica package to the appropriate order where the infinite approx- 
imation leads to the exact solution. 



5 Results and discussion 

Numerical studies have been made for a small amplitude ion- acoustic waves 
in an unmagnetized collisionless plasma consisting of a mixture of a weak 
relativistic warm ion-fiuid and isothermal electrons. We have derived the 
Korteweg-de Vries equation by using the reductive perturbation method 
[13]. The time-FKdV equation is derived from the Eular-Lagrangian us- 
ing Agrawal's method [4]. The Riemann-Liouvulle fractional derivative is 
used to describe the time fractional operator in the FKdV equation. He's 
variational- iteration method [21-23] is used to solve the derived time-FKdV 
equation. 

However, since one of our motivations was to study effects of a relativis- 
tic factor R = uq/cq and time fractional order a on the existence of solitary 
waves. The present system admits a solitary wave solution for any order 
approximation. In Fig.(l), a profile of the bell-shaped solitary pulse is ob- 
tained. Figure (2) shows that both the amplitude and the width of the 
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solitary wave increases with the relativistic factor R. Also, the increasing of 
the time fractional order a decreases the soliton amplitude as shown in Fig 
(3). 

In summery, it has been found that amplitude and width of the ion- 
acoustic waves as well as parametric regime where the solitons can exist is 
sensitive to the relativistic factor R. Moreover, the time fractional order 
a plays a role of higher order perturbation theory in varying the soliton 
amplitude. The application of our model might be particularly interesting in 
some plasma environments, such as ionosphere. 
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Figure Captions 

Fig (1): The electrostatic potential, $(^, r), vs ^ and r, for v — 0.04, 
a = 0.03, a = 0.2 and R = 0.03. 

Fig (2): The electrostatic potential, $(,^,r), vs { at r = 5, v = 0.04, 
a = 0.03 and a = 0.2 for different values of R. 

Fig (3): The amplitude of electrostatic potential, $(0,t), vs q; at = 
0.04, a = 0.03 and R = 0.03 for different values of r. 
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